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Abstract 

The sound mode hydrodynamic dispersion relation is computed up to order for 
a class of gravitational duals which includes both Schwarzschild AdS and Dp-Brane 
metrics. The implications for second order transport coefficients are examined within 
the context of Israel-Stewart theory. These sound mode results are compared with 
previously known results for the shear mode. This comparison allows one to determine 
the third order hydrodynamic contributions to the shear mode for the class of metrics 
considered here. 



1 



1 Introduction 



Hydrodynamics describes the behavior of a fluid on length and time scales which are much 
longer than any microscopic scale. The hydrodynamic stress-energy tensor is constructed 
as a derivative expansion in the fluid velocity, and is required to respect equilibrium ther- 
modynamics and the symmetries in the problem. Such a derivative expansion will contain 
unknown coefficients (transport coefficients). In first order hydrodynamics (an expansion of 
the energy-momentum tensor which contains at most one derivative) without a conserved 
charge, the transport coefficients which enter are the shear viscosity rj, and the bulk viscosity 
(. See [H [2] for a more complete introduction to hydrodynamics. 

The plasma created at the Relativistic Heavy Ion Collider (RHIC) appears to be both 
strongly coupled and well described by hydrodynamics [3 HI E]. Transport coefficients are 
necessary input for hydrodynamic simulations of the RHIC plasma; it is desirable to calculate 
them, but the strong coupling renders conventional perturbative calculations unreliable. 

The AdS/CFT correspondence O [TJ [8] (or more generally, 'gauge /gravity duality' or 
'holography') has become a useful tool in describing strongly coupled systems, and has en- 
joyed arguably its greatest success calculating hydrodynamic transport coefficients. There 
are many ways to calculate the transport coefficients of a strongly coupled gauge the- 
ory using an extra-dimensional gravity dual. One can compute correlation functions of 
the stress-energy tensor and use Kubo formulas or examine the poles of such correlators 
O [ini [HI US [H]. Alternatively, one can examine the behavior of the gravitational back- 
ground under perturbations and determine the dispersion relation for such perturbations by 
applying appropriate boundary conditions. Comparison with the expected dispersion rela- 
tion from perturbations of the energy-momentum tensor yields formulas for the transport 
coefficients [T5l \T6\ [T7]. In addition, the black hole membrane paradigm has been employed 
to calculate the hydrodynamic properties of the stretched horizon of a black hole. In many 
cases, the transport coefficients calculated on the stretched horizon coincide with the trans- 
port coefficients in the dual gauge theory [181 UHl EQl EH [221 ES] . Recently, the work of [23] 
provides yet another way to compute hydrodynamic transport coefficients by deriving the 
equations of fluid dynamics directly from gravity. This work has proved quite influential, 
and has led to much subsequent research [25l [261 [23, 12H1 ISH [30] . We use the gravitational 
perturbation approach similar to [13 [IS], as this work is an extension of a calculation in 



The quintessential example of the success of gauge/gravity duality is the case of the shear 
viscosity. In p3] a formula for rj/s (s is the entropy density) applicable to a wide variety of 
gravitational duals was derived. Application of the formula continually resulted in 



It was later shown that this relation holds for all theories of Einstein gravity [231 [SI] (as- 
suming the dual gauge theory is infinitely strongly coupled). It is quite remarkable that 
this result holds for both conformal and non-conformal theories, and is independent of the 
number of dimensions. In [18], the relation led to the conjecture that rj/s > l/Air for all 
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substances, which has observable consequences for the plasma created at RHIC, even though 
the gravitational dual to quantum chromodynamics (QCD) is not currently known. It is thus 
desirable to find other such universal behavior from gauge/gravity duality in hopes that it 
may have implications for heavy ion collisions at RHIC and at the Large Hadron Collider 
(LHC). 

In the past few years, much work has been done to extend previous analyses to second 
order hydrodynamics [HlESlESlESlEZlESlEniEDlESlEn], which attempts to repair technical 
problems in first order hydrodynamics regarding causality. Most of the work on second 
order hydrodynamics so far has focused on conformal theories. It is notable that a universal 
relation between second order hydrodynamic transport coefficients of a conformal theory was 
presented in [31], though it is not known whether this relation still holds for non-conformal 
theories. 

In this work, we examine second order hydrodynamics for a special class of gravitational 
backgrounds. This class is not necessarily conformal, and includes both Schwarzschild AdS 
(conformal), and Dp-Brane (non-conformal) backgrounds. Specifically, we extend the analy- 
sis of [T7] to the next hydrodynamic order by computing the sound mode dispersion relation 
for such a class of backgrounds. We also discuss implications for second order transport coef- 
ficients within the context of the Israel- Stewart formulation of second order hydrodynamics 
[35]. 

The paper is organized as follows. In section [21 we detail the specific type of gravity dual 
on which we focus in this work. In section [3l we add hydrodynamic perturbations on top of 
this background, and review the calculation of [17], since this work is an extension thereof. 
In section HI we extend the calculation to the next hydrodynamic order. The hydrodynamic 
dispersion relation w{q) in the sound mode can be written 

w{q) = Wiq + W2q'^ + W3q^ + ... (2) 

The main result of this section (and indeed the paper as a whole) is the calculation of w^. 
In section [5l we discuss the implications of our formula for for second order transport 
coefficients within the context of Israel- Stewart theory. In section [6l we compare our results 
to those of a different hydrodynamic mode, the shear mode. It was shown in [33] that 
the shear mode dispersion relation contains contributions from (currently unformulated) 
third order hydrodynamics. Still, by comparing the shear mode and the sound mode we 
can determine the value of these extra contributions for the class of metrics we consider 
here. Finally, in the last section we summarize our results and present prospects for future 
investigation. 

2 Background fields 

Because this work builds upon the paper [17J, we use this section to review the setup and 
notation presented therein. 
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2.1 Black brane background 

In [17] , the sound mode was analyzed for black brane gravitational backgrounds in p + 2 
dimensions. The matter supporting the metric was assumed to be one or more scalar fields. 
The form of the action is assumed to be 



(3) 



S = -r^r- I dF^^x^ [r-IY. d,<t>kd^(l^k - U{(l)^, (t>2..4n)\ • 
lb7rGp+2 J V fc=i / 

The energy-momentum tensor derived from the action is 

n 

MGp+2T^u = ^Yl i(^t^'Pkdu(pk - gi„.^<f>k) (4) 

k=l 

n 

^'i'k = UI^A0fc5^0fe + f^(01,02---0n)- (5) 
k=l 

The metric takes the form 

ds^ = goo{r)dt'^ + gxx{r)dxjdx^ + grr{r)dr^, (6) 

where j = l,2...p, and t is the time coordinate. If this metric is presumed to be dual to a 
strongly coupled gauge theory in p + 1 dimensions, such a theory would live on the boundary 
at r — i> oo. The position of a horizon is assumed at r = rg, and the metric components are 
assumed to behave in the standard way near a black brane horizon, namely that 

c/ooW ~ -7o(r - ro) + 0(r - ro)2 (7) 

g„.{r) ^ ^1^ + 0(1) (8) 
r — tq 

Qxxi^^) ~ 9xx 

(ro) + 0(r-ro). (9) 

The quantities 70, 7r and gxxi^o) are independent of r, though they may depend on tq. The 
Hawking temperature of this metric is 

T = -^& (10) 
Air y 7j. 

For future convenience, we define the function 

F(r) = -g,,{r)g^%r), (11) 
and we use the following notation for the logarithmic derivative 

^L[X{r)] ^ (12) 

The prime denotes derivatives with respect to r unless otherwise noted. Throughout this 
work, our general relativistic conventions are those of [36] . 
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2.2 Restrictions on the metric 



The sound mode gauge invariant equations for the type of backgrounds mentioned above 
were derived in [T3| in much more generahty than is necessary for this work. Instead, we 
will work with a special class of backgrounds; let us assume that the metric is generated by 
a single scalar field, and that the metric components satisfy the following constraints 

grr{r) = c,gUrr'F{r)VL[F{r)]\ (13) 

F(r) = l-( ^^^^] ' . (14) 
V 9xx{r) J 

Here ci, C2 are constant^]. The first of these constraints is a consequence of the fact that the 
particular combination of Ricci tensor components 

Ro -r: = o (15) 

for any metric that is generated by r dependent scalar fields, as shown in [TTj. The second 
of these constraints (HM is imposed to simplify the calculations. It should be noted that the 
Schwarzschild AdS metric and the Dp-Brane metric satisfy both constraints, and the above 
parametrization conveniently allows us to compute the hydrodynamic dispersion relation for 
both of these important special cases. 

The Hawking temperature for this metric is given by the relation 

Cigxx{roY (16) 



2.3 Scalar field and potential 

Before proceeding with the calculation, it is worthwhile to ask what implications the above 
choice of metric has for the scalar field and the scalar potential. The combination of back- 
ground Einstein equations 

g^'Goo - g'^Grr = -SvrGp+a (^7°°Too - g'^'Trr) (17) 

can be simplified to 

</)'(r)2 = 2grr {Rl{r) - K{r)) . (18) 

Explicitly computing the right side for the special metric chosen yields a relationship between 
and g^^. 

0(r) = /€log[fc(r)] +00 (19) 
where 0o is an integration constant, and 



ft: = ±^/|(p+l-2c2). (20) 



^This choice of metric is the same as the example worked out in |17j . with the correspondence C2 — aoj 
and the additional constraint a2 = p ~ clq has been imposed 
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One can determine the form of the potential by considering the following combination of 
background Einstein equations 

^?°°Goo + g"'Grr = -8nG,+2 (^?°°Too + g'-'Trr) = f/(0). (21) 
Explicitly computing the left hand side for our special metric gives 

u{Hr)) = -T^gUr^r'-'gUrr^-"-' = -T^gUr^r'-'gUry-^^'- (22) 

2Ci C2 2Ci C2 

Using (fT9l) . we find 

U{<P) = -T^gUro^' exp f--«:(0 - 0o)l • (23) 
2ciC2 \ p J 

Thus, the sort of metrics we are considering are those generated by a potential which contains 
a single exponential, similar to the Chamblin-Reall backgrounds examined in [3SIES]- 
Note that since the potential must be independent of temperature, it is required that the 
constant C2, and the combination cigxxiroY'^^ must themselves be independent of tq. 

3 Hydrodynamic fluctuations 

One can access the hydrodynamic regime of the dual gauge theory by examining perturba- 
tions of the gravitational background, and following the prescription of [15]. Here we briefly 
review this method 

3.1 Method for determining w{q) 

In the case at hand, one should allow for fluctuations Qf^i, g^u + 6gf^u and 0^0 + 6(t>. The 
space-time dependence of the fields is presumed to be 

5g^,it,z,r) = e^(^^-"*)V(^)> (24) 
S^{t,z,r) = e*("*)50(r). (25) 

Here we use the coordinate z to denote one of the spatial coordinates: z = Xp, and w and q 
are the energy and momentum of the perturbation. For the sound mode, in the gauge where 
h^r = 0, the only non-zero fluctuations are [HI [16]: hoo{r), ^^^Yl^ll hii{r), hzz{r), hoz{r), 
and (50(r). 

Turning on these perturbations, and expanding the background equations of motion to 
first order in the perturbation leads to a set of linearized equations for the perturbations 
mentioned above. To proceed, one should take appropriate combinations of the resulting 
linearized equations and construct equations involving only gauge invariant variables, those 
which do not transform under the diffeomorphism 

V ^ V-Vre.-V(°)^^, (26) 

m -> m-ud'<P) (27) 
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for any vector = ^^(r)e**^'^^~'^*^. (Here, Vlf-* is the covariant derivative with respect to the 
background metric). 

Solving the resulting gauge invariant equations perturbatively in the hydrodynamic regime 
g -C T, and imposing an incoming wave boundary condition at the horizon and a Dirich- 
let boundary condition at the boundary leads to the hydrodynamic dispersion relation w{q). 
One can then compare this dispersion relation to the expected hydrodynamic form to relate 
the transport coefficients to the field components in the gravity dual. 

3.2 First order hydrodynamics for this metric 

The gauge invariant equations which need to be solved for this type of metric were derived 
in [T7j. They involve two gauge invariant variables Zq and Z^. The equations are 



3rr „ 2 172 ' 



-9 



-99 



a 



ip2 



-Z' 



2Z^4)'F adr 



i W 2 



pl^ilgxx] 



-99' 



(28) 



zdr 



where 



a{r) = q (p — 1 



gg"-Z'^] - Z^g„ {q'g'''' + w'g 
VL[goo{r)]\ 



00 



pw 



(29) 



(30) 



VLlgxccir)]] F{r)- 

Solving these equations perturbatively with appropriate boundary conditions leads to the 
dispersion relation w{q). Classically, waves can enter the black hole's horizon, but cannot 
be emitted from there. The standard way to apply this 'incoming wave' boundary condition 
is to make the ansatz [151 [IS] 



Zo{r) 
w{q) 



iw{q) 



F{r)-'i^ (^Fo(r) + qYi{r) + q^Y2{r) + ... 
F(r)-^ (F^o(r) + qY^^{r) + q^Y^^ + ••.) 



Wiq + W2q^ + W'iq^ + 



(31) 

(32) 
(33) 



The functions Y must be regular at the horizon in order for the incoming wave boundary 
condition to be satisfied. One then inserts this ansatz into the gauge invariant equations, ex- 
pands the resulting equation in powers of and solves for the functions Y . Finally, applying 
Dirichlet boundary conditions at the boundary (r oo) gives the dispersion relation. 

First we solve for the functions Y^i. Inserting the above ansatz into (l29l) and expanding 
the result in powers of q leads to 



-99" ( Y;, 



IWi 

2^ 



Vl[F]Y^, 



+ 0{q' 



(34) 
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It should be noted that we have omitted some terms which are proportional to {R^ — R^) 
since they vanish by the background equations of motion. 

Solving the equation for Y^q (and using the fact that \f--gg'"^ oc ^'^[F] from (IT51) ) yields 
a solution 



>0o(r) 



A;o + fci log(F(r)) 



(35) 



but only the constant term can contribute due to the assumption of regularity at the horizon, 
and presumed near horizon behavior of the metric ([7]). Finally the constant must be set to 
zero by the Dirichlet boundary condition at infinity. Proceeding now to higher orders in g, 
one finds that the equations always reduce to the same as that for lo? and as a result = 
to all orders in q. 

Proceeding with the calculation one must now insert the incoming wave ansatz into the 
remaining gauge invariant equation (128|) . 
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J -99 ry, 



0, 



(36) 



solve for the Yi functions, and impose regularity at the horizon and a Dirichlet boundary 
condition at r oo. These steps were completed up to 0{q) in [17]. The results for the 
type of metric we consider here are summarized below. 



Yi{r) 

W2 



( F(r) - 1 ) 

-2(p - C2) 

2ttTp 



(37) 
(38) 

(39) 
(40) 



Here yo and yi are constants. 

Comparing these results to the dispersion relation expected from first order hydrody- 
namics 



Wi 
W2 



V 



p — 1 



(41) 
(42) 



P \ p 

one gains knowledge of Vs (speed of sound), r] (shear viscosity), and C (bulk viscosity). 

rj/s = I/Ati (43) 

(44) 



/2c, 



- 1 



p 



C/v 



(45) 
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Here s = (e + P)/T is the entropy density, e is the equihbrium energy density, and P is 
the equihbrium pressure. Note that the conjectured bulk viscosity bound of Buchel ^0] 
is saturated for metrics of the type we consider. It was also shown in ^/Tj that the above 
results agree with previous calculations for the Schwarzschild AdS metric (with the choice 
C2 = {p + l)/2. These results also agree with the calculation of [16] for the Dp-Brane metric 
with the choice C2 = p{7 — p)/(9 — p). 



4 Solution for ws 

It is the purpose of this work to extend the above calculation to the next hydrodynamic 
order, and thus to determine the next coefficient in the dispersion relation w^. 



4.1 Equation for Y2 

We have already shown that vanishes to all orders in q, so then it remains to return to 
(|36ll . insert the incoming wave ansatz, expand in powers of q, and insert the solutions (|37l 
- Hoi) . Completing these steps, one finds the following differential equation which must be 
solved for Y2: 



dr 



-'2 



F' 



+ 



[l + FY 



Y2 + 



yo 



where 



xq = 2 \ 5wi — 1 



F V 9xx{roF 
2wiW3i'inTy\ 



0. (46) 



(47) 



In writing the above expression, we have replaced all occurrences of the constant C2 which 
appears in the metric with wi due to the relation fl39l) . and have removed the constant ci in 
favor of T using f|T6|) . 



4.2 Solution for Y2 

4.2.1 Solution to the associated homogeneous equation 

The associated homogeneous equation for (H6!) is the same as the homogeneous part of the 
equation for Yi which was considered in [T^. The solution can easily be found by first making 
the ansatz Y2{r) = ?/2a(l — F{r)), and then using the technique of reduction of order once 
this solution is found. The general solution to the homogeneous part of the above equation 
is 



Y2hir) = y2ail - F{r)) + y2b[il - F{r)) log[F(r)] + 4 
Here the subscript h stands for homogeneous, and y2a, y2b are arbitrary constants. 



(4J 
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4.2.2 Particular solution 



In order to find the general solution for I2, we must now find a particular solution to the 
inhomogeneous equation. Since the homogeneous solution is known, one can construct a 
particular solution using the method of variation of parameters. For completeness, this 
method is outlined in Appendix \M Applying the method to the case at hand gives the 
following solution (details are given in Appendix [Bl) . 



Q{r) 



yo 



,(F(r) 



fl-F{r)Y F'(f) 
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1 + F{r)\ F'{r] 



1 - F{r) J F{r) 



Q{r) dr 



wf + 



9xx{r) 



F{r) - wl 



dr. 



l + F{r)J F{r) [ ' ' fc(ro)f 
Thus, the general solution for Y2 is given by 

Y2{r) = Y2h{r) 

To proceed, it is convenient to change coordinates. We define the coordinate u by 



u 



gxx[ro) 
, 9xxir) 



gxx[ro) 
, 9xxir) 



(49) 
(50) 

(51) 

(52) 
(53) 



so that 

F{u) = \-u^, 

and the horizon is now located at m = 1, and the boundary is located at m = 00 In terms of 
the new coordinates, Y2 is given by 



Y2M 



y2aU + y2b 

yo 



n^log[l - n^] + 4 



{2-u 



2\2 



uHi 



W 



2+ ^-4/(1+^,2) 



(2-t;2)2(l -t;2) 



(54) 

dv du 
(55) 



The inner integral can be written in terms of Hypergeometric functions, but the result is not 
particularly enlightening, so we do not reproduce it here. The full analytical form of this 
function is not needed to determine the dispersion relation. 



4.3 Boundary conditions 
4.3.1 Regularity at horizon 

The first boundary condition which must be applied on Y2 is the condition of regularity at 
the horizon. In order to do so, one must extract the coefficient of the logarithmic divergence 

^Here the assumption is that C2 > 0, and that gxx{f — * 00) ^ r" with n > These assumptions hold for 
the Schwarzschild AdS metric for any positive p, and for the Dp-Brane metric provided p < 7. 
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in the particular solution Y2p. To do so, we first expand the integrand in powers of {u — 1), 
and look for the coefficient of the {u — term. After integration, this term will lead to 
the logarithmic divergence. 

With the aid of Mathematica, we find the nested integral can be expanded near the 
horizon as 



/"- 



w 



(2-i;2)2(l-t;2) 



dv 



1 ( 1 - 3w;f 2 

2 T^^-"-? 



2 + ivr — Hn 



-2 



l + wi 



+ 0{u-l) 



where Hn{a) is the 'Harmonic Number' defined as 



HJa 



dx. 



X 



Using this expansion in I2, we find that near the horizon, 

Y2p{u ^ 1) 



du 
4(1 --u) 



1 - 3u;2 . ( . 

^ - wl [2 + t7r-Hn 

l-wf V 



Going back to the general solution for I2, one thus finds 
Ri log[l - u] \ y2b + 



yo 



1-wi 



The requirement of regularity ai u = 1 thus gives 

yo 



1 - 3wf 2 L ■ 



i + wi 



y2b = 



wi\2 + i7r - Hn 



l-3wf 



(47rT)2 

For future convenience, we make use of the identity 
to write 



2 



1 



1 + a 



Hr, 



l-2wf- 3wf 



1 + wlJ 2wl{l-wl) 
Substituting this into the equation for we find 



y2b 



yo 



wl[2 + iTr- H„ 



(56) 



(57) 



(58) 



+ 0{1) 



(59) 



+ 0{1). 



(60) 



(61) 



(62) 



(63) 
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4.3.2 Dirichlet boundary condition at m = 



Finally, we must apply the Dirichlet boundary condition at m = 0. We proceed as above, by 
expanding the integrand of near u = 0. Mathematica gives 



u V 



w 



2 + y-4/(l+«,?) (1 



V — wt 



yo 



(2 -t;2)2(l 

so that 

Y,(u ^ 0) ^ {Ay2t + Oi 

Doing the integral, one finally has 

Y^iu 0) ^ 4y2b - 

and applying Y2{u ^ 0) = gives 

42/2fe 



dv 



w 



i(2 + z7r) + 0(M), 



Xq — 8u Wi{2 + in) 



+ 0{u 



-2^ 



du 



yo 



xo + Awl (2 + «vr) + 0{u) 



yo 



Using dM]), 



Xq + Awf (2 + in) 



2wf 



Xo = -4 

Finally, fH7j) allows us to solve for w^- The result is 

wi{l-wj) 



W3 



(47rT)^ 



(64) 



(65) 



(66) 



(67) 



(68) 



(69) 



This is our main result. To summarize, we have computed the coefficient of the term in 
the sound mode hydrodynamic dispersion relation for a specific class of metrics (see f ll4p ). 
This class of metrics contains two constants which we denote ci and C2 . Our expression for 
1^3 should necessarily depend on these constants, but we have eliminated ci in favor of T 
using f|T6|) and C2 in favor of wf = due to the relation 



5 Transport coefficients in Israel- Stewart theory 

Now that we have computed the dispersion relation to 0{q^) in the sound mode, we are in a 
position to examine the implications for second order hydrodynamic transport coefficients. 
Second order hydrodynamics attempts to repair some technical problems regarding causality 
in the first order theory. This subject was first broached by Miiller [H], and later by Israel 
and Stewart [35]. Recently, the formulation of second order hydrodynamics presented in 
[33] has gained popularity, though at present it is only applicable to conformal theories. (It 
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should be noted that recently, some progress has been made in generalizing the work of [33] 
to non-conformal theories [12] )• The metrics we consider are not necessarily conformal, and 
thus we will use the Israel-Stewart formulation. 

Israel introduced five new transport coefficients that appear in the hydrodynamic ex- 
pansion of the energy momentum tensor. In what follows, we use the same notations and 
conventions as p52j. Three of these five transport coefficients are relaxation times associ- 
ated with the diffusive, shear, and sound mode, and are denoted by (rj, r^, rn) respectively. 
There are two other transport coefficients which are related to coupling between the different 
modes, ao,ai. 

In [32], the sound mode dispersion relation was computed in terms of these transport 
coefficients. 



Wi 

W3 



i 

r7 



V 



2v.Ts 



p — 1 
P 



p J Ts 



+ 



i) 




- ^ )1 


pj 


Ts 


ATsj _ 



(70) 
(71) 

(72) 



These relations were derived within the context of a 'decoupled ansatz' which presumes the 
background contains no R-charge. The cases which we will consider below (Schwarzschild 
AdS and Dp-Brane) fit this criteria. Our backgrounds are generated by scalar fields only, 
and thus any gauge field necessary to provide R-charge is absent from the cases we consider 
here. A more complete list of assumptions regarding this dispersion relation can be found 
in t32j. 

Comparing fl72|) to our main result (l69l) . and eliminating 77 and C, from the relations 
fl43|45l) gives the relation 



(p-i) 



P 



i8TrT)vUp-l] 



l + vt 



l + 2H„ 



2vl 







(73) 



As expected, the coefficients and rn cannot in general be determined separately using this 
method. Still, if one of these coefficients is known, the above relation allows us to determine 
the other. 

Let us now explicitly check that our results agree with other calculations in the case of a 
conformal background. The Schwarzschild AdS black hole metric in the near horizon limit 
takes the form 



Fir) 



r 

1? 



F{r)dt^ 

)' 



+ 



r'^F{r) 



(74) 
(75) 
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where L is the radius of curvature of the AdS space. For this metric, the parameter C2 
(p + l)/2 and thus = 1/p. In this case, ( 1751) gives 



SAdS 



4nT 



p+1 



1 



(76) 



which is in ag reement with [SB], p7|2|. This is a non-trivial check on our calculation; the 
cited results were arrived at by completely different methods than those we employ here. 
Furthermore, it should be noted that (1761) confirms a conjecture made by Natsuume in \22fi 
Finally, we can also consider the case of the Dp-Branes. In the Einstein frame, the metric 
can be reduced to |18l |16| 



F{r) 



9-p 

V 



-F{r)dt^ + dxjdx^ 

7-p 



+ 



p^-8p+9 



dr"^ 



(77) 
(78) 



for this metric, the parameter C2 = (7p — p^)/(9 — p), and v1 = (5 — p)/(9 — p) \W\. Inserting 
this into (1731) one finds 



(9-p)(l-p)r,^^ = (3-p)V^ 



TcT (5 — p) 



P 



7 — p \7 — p, 



(79) 



This formula agrees with previous computations for p = 1 and p = 4 in 



6 Comparison with the shear mode 

Though we have exclusively dealt with the sound mode in this work, we now compare our 
results to similar ones from the shear mode. Let us define the shear mode dispersion relation 
in the same way as [32|, 

U'(g)shear = 'iD^q^ - ^D^T^hearg' + O(g') (80) 

where 

^« = — = (81) 

A formula for Tghcar was given in pT] which is applicable to a wide variety of metrics, including 
the special metrics we have considered in this note. It states that 



^(^o) 



''"shear i / dv 



9 



rr \ 



^ ( Dir) 



(82) 



■^In comparing with the resuhs of ^21.1 , one needs to employ the identity ([6T1) to see the agreement 

''To be precise, the conjecture is confirmed for the case of p > 2; the case of p = 1 should probably be 

checked separately as the derivation of the gauge invariant equations in [TB], [T7] rely on at least 2 spatial 

dimensions. See |43| . where first order hydrodynamics is examined for p ~ 1 
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where 



Dir) . ^ r dr'^pk^^ . (83) 

V-^oo(r)^rr(r) ^-9{r')gxx{r') 



Using the special metric fll3tll4p . and the relationship (fT6|) in this formula yields 

r,.„. ^ (2 - I) ^ (^) . (84) 

This result agrees with the special cases of the Dp-Brane and the Schwarzschild AdS metrics 
as computed in [22] . 

Previously, it was thought that Tghcar = Tw, but recently the authors of [33J showed 
that this is not the case because the coefficient of the term contains not only r^, but 
also contributions from (currently unformulated) third order hydrodynamics. We can now 
determine these unknown contributions for metrics which obey (I13fl4l) . 

Let us parametrize 

Tshear = + A (85) 

where A denotes the unknown contributions from third order hydrodynamics. Combining 
( IH^ and (!73|) allows one to solve for A. Evidently, 

The correction A does not appear to be universal in the sense that the first term is not 
present in the case of a conformal theory. In the conformal case, A is still not 'universal', 
because it depends on the number of dimensions of the theory. In the future, when the 
particular transport coefficients which comprise A are known, it will be interesting to see 
whether there is any universal relationship between these unknown coefficients, r^, and tjj- 
We can easily check that the formula fl5^ reproduces the results in the well known SAdS^ 
metric, which is dual to AT = 4 supersymmetric Yang-Mills theory at finite temperature. In 
this case, p = 3 and v"^ = 1/3. Immediately, we have 

^SADS5 ^ L. (87) 



In [321 133], it was found thal| 



27rT 



SADS5 ^ 1 - log(2) 
' shear 2,7tT 

.SADS5 _ 2-log(2) 



^.Au., = (89) 

It is clear that the relaxation time computed from the shear mode and the sound mode differ 
by the amount predicted by the formula (|86|) ; our results are in agreement with [32l [33] . 



^Even though the formalism of Baier et al. [32] is different than Israel-Stewart, one can check that the 
sound mode dispersion relations coincide in the limit of conformal theories (C — > 0). In an unfortunate clash 
of notations, the relaxation time introduced by Baier et al. is denoted by rn, its Israel-Stewart counterpart 
is T^. 
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7 Conclusion and outlook 



In this work, we have extended the calculation of [171 the next hydrodynamic order in q. 
The main result (169|) is the coefficient of the term in the sound mode dispersion relation. 
This result is applicable to metrics which obey (fT3] - [T^ . These metrics are not necessarily 
conformal, and contain an arbitrary number of spatial dimensions p with p > 2. 

Information about second order transport coefficients was presented within the context of 
the Israel-Stewart theory (specifically within the formulation presented in [32]). In general, 
a relationship fl73l) between two transport coefficients r^r and rn can be determined. 

In the conformal case of the Schwarzschild AdS metric, the relation mentioned above 
allows the determination of the coefficient r^. We have verified that our results agree with 
those calculated from different methods. 

Finally, by comparing the sound mode dispersion relation to the shear mode discussed 
in [211 122], we were able to determine the contribution of third order hydrodynamics to the 
shear mode (it was pointed out that such contributions would be present in [33]). 

As mentioned in the introduction, it is desirable to find other universal relations among 
transport coefficients, as such relations sometimes lead to observable consequences. The 
main results of this paper (l69l) . (1731) and the third order hydrodynamic contributions to the 
shear mode fl86|) are applicable to certain gravity duals which may or may not be conformal. 
However, despite the fact that these relations appear to be applicable to many theories, they 
do not seem to be universal in the same way a.s rj/s = l/47r, and the relation presented in 
[51] . For example, the number of spatial dimensions p enters explicitly into our formulas, 
whereas a universal relation should not depend on this quantity. 

Furthermore, we stress that the class of theories examined herein is quite limited. It 
seems likely that any generalization will explicitly contain the bulk viscosity (, which we 
were able to eliminate due to the relation fH5|) . 

Of course, it would be interesting to re-examine the results here for a broader class of 
gravity duals. Unfortunately, the gauge invariant equations derived in are difficult to 
solve analytically except in the special case presented here, though perhaps insight could 
be gained by approaching the problem numerically. It would also be useful to attempt to 
generalize the results of |T7] by including other matter fields in addition to the scalar fields 
considered there. It seems unlikely (though perhaps not impossible) that similar special 
cases would be analytically solvable after the addition of different kinds of matter. These 
are issues which should be investigated in the future. 
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Appendices 

A The method of variation of parameters 

One fundamental technique used in generating particular solutions to inhomogeneous dif- 
ferential equations is the method of 'Variation of Parameters'. The method can be used to 
find a general solution to a second order linear differential equation if the solution to the 
associated homogeneous equation is known. 

The theory behind the method can be found in any textbook on differential equations. 
Here, we simply present the essential formulas. Consider a differential equation 



and assume that the functions yi{x) and y2{x) are linearly independent, and satisfy the 
associated homogeneous equation. That is, 



and similarly for y2{x). It can be shown that the function yp{x) is a solution to the inhomo- 
geneous equation, where 



y"{x) + p{x)y'{x) + q{x)y{x) = g{x) 



(90) 




(91) 




(92) 



Here, W{x) is the Wronskian 




(93) 



To arrive at the form of the function yp used in the text (l50i) . let us define 




(94) 



Then, 



W{x) = h'{x)yl{x) 



(95) 



and 




(96) 



One can see that this is equivalent to 




(97) 



by performing an integration by parts. 
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B Particular solution for Y2{r) 



One can now apply the method detailed in Appendix |A] to the differential equation for 1^2- 
Using the notation from Appendix |Al (and changing independent variables from a; to r), we 
have 



h{r) 
9{r) 



l-F(r), 
log[F(r)] 



1 -F(r)' 



/\2 



yo 



F{l + F) {AnTy 



( 2 ^ gUry 

Xo + — — + 



After a bit of work, one can find the particular solution 



y2p{r) 
Q{r) 



Vo 



[I -Fir)) 



'1 + F(r)\ F'ir] 



1 - F{r) J F{r) 



:i+F{z)r 



Xo + 



1 - F{z) 
F{z) 



gxxiroY 



Q{r) dr 



gxxiro)P 



(F-wI 



(98) 
(99) 

(100) 



F{z)-wl 



(101) 

dz. 
(102) 



The term involving xq can be integrated directly, by changing variables from z to F. This 
simplifies the solution to 



Q{r) 



+ 



yo 



(47rT)2 (AnT) 
(I- F{z)V F'{z) 



-XF{r) - 1) 



1 + F(r)\ F'(r) 



l-F(r); F(r) 



Q{r) dr 



\l + F{z)J F{z) 
which is the form presented in the text. 



, gxxKZf („, . 2\ 



gxxiro) 



dz, 



(103) 
(104) 
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